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Abstract 

This work gives a van’t Hoff law expression of Langmuir constants of different species for deter¬ 
mining their occupancy in the nanocavities of clathrate hydrates. The van’t Hoff law’s parameters 
are derived from a fit with Langmuir constants calculated using a pairwise site-site interaction po¬ 
tential to model the anisotropic potential environment in the cavities, as a function of temperature. 
The parameters can be used for calculating clathrates compositions. Results are given for nineteen 
gas species trapped in the small and large cavities of structure types I and II [Tt. The accuracy 
of this approach is based on a comparison with available experimental data for ethane and cyclo¬ 
propane clathrate hydrates. The numerical method applied in this work, was recently validated from 
a comparison with the spherical cell method based on analytical considerations [1 


1 Introduction 

A clathrate is an ice-like crystalline solid consisting of water molecules forming a cage structure around 
smaller guest molecules under suitable conditions of low temperature and high pressure. On Earth, it 
is considered that clathrate hydrates are the most important reservoirs of fossil energy [Hi, and that 
favourable conditions for gas hydrate formation exist in about 25% of the earth’s land mass. Moreover, the 
thermodynamics conditions of pressure and temperature prevailing in the oceans are such that hydrates 
should easily be formed in about 90% of the ocean or sediments. The most common guest molecules in 
terrestrial clathrates are of organic aliphatic nature like methane, ethane, propane or butane, but other 
small inorganic molecules like nitrogen, carbon dioxide, and hydrogen sulfide can also be trapped in the 
cages of clathrates [2 El El 0 IS M HD]- In the advent of global warming, these clathrates can enhance 
the temperature rise when the trapped species are released. Clathrate hydrates are also suspected to 
be extensively present on several planets, satellites and comets of the Solar System. Planetologists are 
thus concerned with the possible clathrate impact on the distribution of the planet’s volatiles and on the 
modification of their atmosphere’s compositions mi. Hence, it is of great interest to correctly determine 
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the amount of species potentially trapped in the cages of clathrates, i.e. the fractional occupancy of guest 
species under the thermodynamic conditions (pressure and temperature) prevailing in the regions where 
clathrates might form. 

From a theoretical point of view, the thermodynamics of the formation or dissociation of clathrates 
is most often based on the model developed by van der Waals and Platteeuw [12] following the same 
hypotheses under which was developed the adsorption theory of Langmuir [T3]. The Langmuir isotherms 
of adsorbed molecules on a surface are determined from the calculation of the Langmuir constant, which 
is also the main parameter to be considered in the determination of the amount of species trapped in the 
clathrate cages as a function of pressure and temperature. 

To calculate these Langmuir constants, most of the models are based on a molecular description of 
the guest-water interactions using a Lennard-Jones or Kihara potential form. The parameters of these 
potentials are usually empirically obtained from experimental data of phase equilibrium. Such models 
most often neglect interactions of the guest molecules with water beyond a few cages only and are therefore 

questionable na Hi m ua mm • 

Moreover, it is generally assumed that the environment of the cage in which a gas molecule is trapped 
in clathrates is of spherical symmetry. Whereas this assumption may be justified for molecules such as 
CH 4 or NH 3 , it is certainly not well-suited for molecules such as CO 2 and N 2 O or SO 2 which are of 
cylindrical or oblate symmetry and for which the hypothesis of a free rotation at the center of a spherical 
cage is no longer valid. 

In the case of ethane and cyclo-propane clathrate hydrates, an analytical method based on the spher¬ 
ical cell model has been used to extract spherically averaged intermolecular potentials from experimental 
data on the temperature dependence of the Langmuir constant by Bazant and Trout ]20] . For other guests 
in clathrate hydrates, in particular, argon, hydrogen, nitrogen, methane, ethane, propane, cyclo-propane, 
and carbon dioxide, other workers mmmm have explicitly taken into account the angle-dependence 
of the guest-water intermolecular potential in an atom-atom or site-site description to calculate the 
corresponding Langmuir constants. However, in most of these studies the water and guest molecules 
were simply described as one, two or three interaction sites for the Lennard-Jones or Kihara potential 
contributions although more than three sites are involved. 

This work aims at providing a van’t Hoff law expression of the Langmuir constant for single guest 
molecules incorporated in clathrate hydrates as a function of the temperature by improving the potential 
model used in the determination of the Langmuir constants , that is by using an atom-atom and site- 
site potential and by considering explicitly the effect of water molecules beyond the trapping cage and 
the resulting anisotropic environment. In the present work, the Langmuir constant is determined by 
taking into account all the external degrees of freedom of the guest molecules, i.e. the center of mass 
(c.m.) translational motion and the orientational motion in a true crystallographic clathrate lattice, not 
necessarily of spherical symmetry as it is often assumed when using the van der Waals and Platteeuw 
model P2]. In the following, we recall in section 2 the model used for the calculations of the Langmuir 
constants that are necessary to determine the fractional occupancy of guest species in clathrates. Then, 
in section 3, the geometry and interaction potential considered here are described. Finally, in section 
4, the coefficients for a simple van’t Hoff expression of the Langmuir constants are given for a large set 
of guest molecules. The Langmuir constants calculated using this model are compared with available 
experimental data, i.e., for ethane, cyclo-propane guest molecules. 

2 Statistical thermodynamic approach 

In contrast to natural ice which solidifies in the hexagonal structure, clathrate hydrates form, as water 
crystallizes, in the cubic system in several different structures which are characterized by specific cages of 
different sizes. The two most common types are ’’structure I” and ’’structure II”. In structure I, the unit 
cell is made of 46 water molecules forming 2 small (12 pentagonal faces 5 12 ) and 6 large (12 pentagonal 
and 2 hexagonal faces 5 12 6 2 ) cages, while in structure II, the unit cell is made up of 136 water molecules 
forming 16 small (12 pentagonal faces 5 12 ) and 8 large (12 pentagonal and 4 hexagonal faces 5 12 6 4 ) cages 
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Calculations of the relative abundances of guest species incorporated in a clathrate lattice structure of 
type I (si) or II (sll) at given temperature-pressure conditions can be performed using classical statistical 
mechanics which allows the macroscopic thermodynamic properties of the clathrates to be determined 
from the interaction energies between the guest species and the clathrate water molecules. 

In 1959 van der Waals and Platteeuw [T2J developed a model of clathrate’s formation in which the 
trapping of guest molecules in nano-cages was considered to be a generalized case of the three-dimensional 
ideal localized adsorption. 

Their model is based on the following hypotheses: 

1. The contribution of the host molecules to the free energy is independent of the occupational mode 
in the cages. This implies in particular that the guest species do not distort the trapping cage. 

2. The encaged molecules are localized in the cages, each of which can never hold more than one guest. 

3. The mutual interaction of the guest molecules is neglected, i.e., the partition function for the motion 
of a guest molecule in its cage is independent of the other guests. 

4. Classical statistics is valid, i.e., quantum effects are negligible. 

From the configuration partition function and the thermodynamic equilibrium condition on the chem¬ 
ical potentials of the guest and host molecules in coexisting phases in clathrate [2T| , the fractional occu¬ 
pancy of a guest molecule K in a given ’’structure type-cage size” t (t = structure-type I or Il-small or 
large cage) can be written as: 


VK,t ~ i + E 


(i) 


where the sum in the denominator includes all the species present in the initial gas phase, Cx,t is the 
Langmuir constant of species K in the structure type-cage size t, and fx is the fugacity of the species K 
which depends on the total pressure P of the initial gas phase and on the temperature T. 

The Langmuir constant in equation (1) depends on the temperature T and on the strength of the 
interaction energy between the guest species I\ and the water molecules in the cage. It is expressed as: 


CK,t = 



VkA r,») 

k B T 


drdCl. 


( 2 ) 


In this equation V^.t (r, Li) is the interaction potential energy experienced by the guest molecule for 
a given position vector r of its center of mass with respect to the cage center and its orientational vector 
$7, and k B is the Boltzmann constant. The integral value must be calculated for all external degrees of 
freedom of the guest molecule inside the structure type-cage size t. 

To evaluate the Langmuir constant Cx,t, two additional assumptions are often made, namely: ( i ) the 
symmetry of the guest molecule’s environment is considered to be spherical and (ii) the guest molecule can 
freely rotate in the corresponding spherical cage (spherical cell potential approximation), in accordance 
with the Lennard-Jones and Devonshire [22) theory applied to liquids. 

Then, the Langmuir constant can be cast as: 


Cx.t = 




(3) 


where R c is the radius of the spherical cage and V B ' t (r) is the spherically averaged potential energy 
between the guest molecule and the clathrate water molecules. 

Note that Eq. (3) is commonly used by planetologists and as a result, it may introduce significant 
inaccuracies in the evaluation of the relative abundances of gas species in clathrate hydrates when the 
assumed spherical symmetry of the interaction potential is questionable. 
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In contrast, the Langmuir constant for enclathrated single guest molecules is determined here by 
taking into account all the external degrees of freedom of the enclathrated molecules, i.e. the center of 
mass (c.m.) translational motion and the orientational motion in a real crystallographic clathrate lattice, 
not necessarily of spherical symmetry. As a consequence, Eq. (3) is no longer valid in our approach and 
should be replaced by a more general form, as explained below. 


3 Potential energy model 

Let us consider a non-vibrating single gas species (atom or molecule) trapped in a small or large cage 
of a clathrate structure of type si or sll. The positions of the hydrogen and oxygen atoms are given by 
their coordinates in a reference frame fixed to the clathrate lattice, while those of the guest molecule 
are usually given in a frame tied to the molecule with the z axis coinciding with the rotational axis of 
highest symmetry as shown for the trapping of a guest molecule (CH 4 for instance) in the small cage of 
the clathrate structure type I, in Figure 1. 



Figure 1: Trapping geometry of a guest molecule (CH 4 for instance) in a small cage of the structure type 
I. (0,X,Y,Z) represents the absolute frame tied to the clathrate lattice. 

The interaction potential energy Vx,t between the guest molecule K and the surrounding water 
molecules of the structure type-cage size f, considered as a rigid clathrate crystal, is modeled as a sum 
of a 12-6 Lennard-Jones (LJ) pairwise atom-atom potential characterizing the repulsion-dispersion con- 
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tributions and an electrostatic part due to charge-charge interactions between the charges in the guest 
molecule and those in the water molecules of the clathrate system. It is expressed as: 


N C Nw Nk 

^ = £££ 

k—1j= 1 i—1 



i qiqj 
47re 0 kyj 


(4) 


where r ij k is the distance vector between the ?th site, with electric charge qi, of the guest molecule (Nk 
sites) and the jth site, with electric charge qj , of the fcth water molecule (Nyy sites) of the clathrate 
matrix (containing Nq water molecules); e t j and are the mixed LJ potential parameters, using the 
Lorentz-Berthelot combination rules as: 


and 


+ a j 


(5) 


where and <r,; are the Lennard-Jones parameters for the i — i interacting atomic pair, which are taken 
from the literature [23 HU US]. The effective electric charges for the water molecules q B = + 0.4238 e, qo 
= - 0.8476 e are from reference [25] and those for the guest species are from references [2511251 i2Tj [251 [251 . 


Table 1: Pure Lennard-Jones potential parameters^ used in our calculations. 



H - H 

0 - 0 

C - C 

N - N 

S - S 

Ne - Ne 

Ar - Ar 

Kr - Kr 

Xe - Xe 

e, (K) 

8.59 

57.41 

42.88 

37.29 

73.79 

43.16 

126.90 

179.85 

226.03 

<Ti (A) 

2.810 

3.030 

3.210 

3.310 

3.390 

2.764 

3.405 

3.650 

3.970 


(a) Obtained from Refs. [221 [2U [25] . 


Calculations of the Langmuir constant Cx.t from Eq. (2) require an explicit determination of the 
external degrees of freedom of the guest molecule, that is, its center of mass (c.m.) position vector r and 
its orientation vector f2 = (tp,9,x). Therefore, we define an absolute frame (0,X,Y,Z) connected to the 
clathrate matrix. 

Figure 2 gives the geometrical characteristics of the guest-water interacting molecules, and a descrip¬ 
tion of the internal positions of the sites in the guest molecule with respect to its frame (G,x,y,z) and its 
external orientational and translational degrees of freedom with respect to the absolute frame (0,X,Y,Z). 

Then, the distance vector r, (;fc in Eq. (4) can be expressed in terms of the position vector r of the 
c.m. of the guest molecule and r, Jk of the jth site of the fcth water molecule with respect to the absolute 
frame (0,X,Y,Z), and of the position vector r/ ? ; of the Ah site of the guest molecule with respect to its 
associated frame (G,x,y,z) (see Figure 2), as: 


r m = r A -r-m. (6) 

It should be noted that the explicit dependence upon the angular degrees of freedom ip, 9 , and x of 
the guest molecule requires the determination of the site position vectors { 77 } with respect to the absolute 
frame. 

Then, by assuming (i) that the guest species and water molecules are rigid, and (ii) that the clathrate 
lattice is undistorted (static lattice), the Langmuir constant is calculated from a six-dimensional config¬ 
urational integral written as: 


Ck.i = 



VK,t(x,y,z,ip,9,x) 

k B T 


dxdydzdip sin 9d9dx- 


(7) 
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Figure 2: Geometrical characteristics of a guest molecule (CH 4 for instance) interacting with a water 
molecule of the clathrate matrix. (0,X,Y,Z) and (G,x,y,z) represent the absolute frame tied to the 
clathrate matrix and the molecular frame, respectively. 


Of course, it should be noted that ( i ) for atomic species, there are no orientational variables, and (ii) 
for linear molecules there is no spinning variable 


4 Results and discussion 

As already described above, we herein use a pairwise atom-atom Lennard-Jones and a site-site electrostatic 
potentials for calculating the interaction of single guest molecules with clathrate matrices containing 4 x 
4x4 unit cells (up to 2944 water molecules) for the structure si, and 3x3x3 unit cells (up to 3672 
water molecules) for the structure sll. Indeed, it has been previously shown that these number of water 
molecules are large enough to ensure a good convergence of the corresponding calculations [501 . Moreover, 
in order to determine the Langmuir constant and its temperature dependence for a guest molecule K 
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in a structure type-cage size t of clathrate, the potential energy surface V/f,t(a;, y, z, ip, 9, \) in Eq. (7) 
is first numerically computed for all the center of mass positions and orientations of the molecule with 
respect to the absolute frame as defined, and then the numerical integration is performed. Note that the 
minimum of the potential energy and the associated molecular position and orientation can also be 
determined in this way. 

4.1 Calculation and fit of the Langmuir constants 

A clathrate is characterized by its degree of occupancy i.e. the fraction of cage sites that is occupied by a 
guest molecule, which can be calculated from the knowledge of the Langmuir constant Ckj- Because such 
constants are not necessarily extracted from experiments in a large range of temperatures, we propose 
here to express them as a function of temperature through a simple van’t Hoff law : 

C K ,t{T) = A K>t exp (B K>t /T), (8) 

in which A (Pa -1 ) and B (K) are constant values fitted to values derived from an atomic description 
of the clathrate systems. This simple van’t Hoff expression could thus be readily used, for instance, by 
planetologists in the large range of temperatures encountered in the Solar System. 

In this work Langmuir constants have been calculated for eighteen species, in the temperature range 
between 50 K and 300 K, as given in Tables 2, 3 and 4 for rare gas atoms, linear molecules and non-linear 
molecules, respectively. 

Table 2: Parameters Ax.t (Pa -1 ) and Bk-,t (K) for the van’t Hoff expression of the Langmuir constant 
for simple guest clathrate hydrates. 


Structure type 

si 

si 

sll 

sll 

Cage size 

small cage 

large cage 

small cage 

large cage 

Guest species K 

Ax,t 

Ax,t 

Ak,i 

An,t 


Bx,t 


Bx,t 

Bx,t 

Ne 

1.0308xl0 -9 

2.7194xl0 -9 

0.8154xl0 -9 

5.6151x 10 -9 


1187.948 

1015.862 

1233.898 

898.062 

Ar 

1.5210xl0 -1 ° 

7.7829xlO -10 

1.0456xl0 -1 ° 

2.4531xl0 -1 ° 


2961.545 

2521.758 

2977.025 

2195.964 

Kr 

0.5985xl0 -1 ° 

3.8799xlO -10 

0.4004xl0 -1 ° 

16.2508xl0 -1 ° 


3885.383 

3454.028 

3789.957 

3021.690 

Xe 

1.9389xl0 -11 

13.8509xl0 -11 

1.2913xl0 -11 

83.2357xl0 -11 


4547.654 

4572.732 

4085.506 

4103.672 
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Table 3: Parameters Ax t t (Pa 1 ) and B Kt (K) for the van’t Hoff expression of the Langmuir constant 
for simple guest clathrate hydrates. 


Structure type 

si 

si 

sll 

sll 

Cage size 

small cage 

large cage 

small cage 

large cage 

Guest species K 

Ax,t 

Ax,t 

Ax,t 

Ax,t 


Bx,t 

Bx,t 

Bx,t 

Bx,t 

h 2 

4.7301 xlO" 9 

16.0695xl0" 9 

5.5295xl0" 9 

64.0074xl0" 9 


1265.757 

1515.721 

1203.620 

873.259 

o 2 

1.3153xl0" 9 

7.4606xl0" 9 

0.8004xl0" 9 

25.8421xl0" 9 


2917.693 

2558.746 

3044.536 

2238.052 

n 2 

3.9496xlO~ 10 

25.6897xl0 _1 ° 

4.8836xlO" 10 

201.3238xl0" 10 


2869.400 

2680.372 

2679.423 

2226.480 

CO 

2.5937xlO~ 10 

16.5331 xlO" 10 

4.5198xl0 _1 ° 

229.8631 xlO" 10 


3518.021 

3075.059 

3088.930 

2275.803 

co 2 

7.7765xl0~ 12 

520.5579xl0 -12 

7.9970xl0" 12 

6907.0012xl0~ 12 


2976.629 

4674.690 

2277.757 

3370.363 

HCN 

7.6653xl0~ 12 

131.1607xl0" 12 

13.9141 xlO" 12 

8224.5133xl0~ 12 


4085.369 

4328.556 

2593.031 

2640.868 

c 2 h 2 

0.9702xl0~ 12 

221.9622x 10 -12 

0.2439xl0 -12 

3909.8421 xl0~ 12 


735.205 

3076.356 

321.114 

2837.467 


4.2 Comparison with available experimental Langmuir constants 

4.2.1 Ethane and cyclo-propane simple hydrates 

In Tables 5 and 6, the Langmuir constants calculated here are compared to experimental values obtained 
from the dissociation pressure data by Sparks and Tester [1] for temperatures ranging between 200 K and 
290 K and between 240 K and 290 K, for ethane and cyclo-propane trapped in large cages of structure si, 
respectively. Moreover, for comparison the Langmuir constants calculated from a simple spherical model 
as often used in the literature are also reported in Tables 5 and 6. 

It can be seen that the values calculated here are in good agreement with those extracted from 
experimental data. The ratios C ca i/C eyi p vary between 0.5 and 1.2 for ethane, and between 0.7 and 1.1 
for cyclo-propane, depending on temperature. Note that the same ratios calculated within the same 
temperature range for the large cages of structure sll lead to values ranging from 320 to 540 for ethane, 
and from 920 to 1100 for cyclo-propane. These results are consistent with the fact that ethane and cyclo¬ 
propane form si-large cage clathrate hydrates only. Moreover, Tables 5 and 6 show that the best agreement 
between the values calculated here and those extracted from experiments is found when temperature is 
higher than 260 K. Below this temperature, the calculated values are systematically lower than the 













Table 4: Parameters Ax t t (Pa 1 ) and B Kt (K) for the van’t Hoff expression of the Langmuir constant 
for simple guest clathrate hydrates. 


Structure type 

si 

si 

sll 

sll 

Cage size 

small cage 

large cage 

small cage 

large cage 

Guest species K 

Akj. 

Ax,t 

Ax,t 

Ax,t 


Bx,t 

Bx,t 

Bx,t 

Bx,t 

h 2 s 

2.3444x 10 _1 ° 

7.2080xl0 _1 ° 

3.6415xlO" 10 

758.3575xl0 _1 ° 


4463.910 

4073.045 

3073.324 

2495.937 

S0 2 

1.2311xl0~ 12 

75.0641 xlO -12 

6.1515xl0" 12 

17926.7530x 10~ 12 


4374.084 

6272.810 

1548.504 

4139.948 

nh 3 

8.6697X10- 11 

27.5220xl0' n 

38.7320x 10” 11 

5031.4000xl0 -11 


5197.361 

4975.753 

3334.132 

2484.181 

ch 4 

8.3453xlO -10 

116.6313xl0 _1 ° 

5.4792 xlO" 10 

829.8039xl0 _1 ° 


2901.747 

2959.901 

2546.660 

2629.194 

C 2 Hg a) 

_ 

3.5164xl0 -11 

_ 

727.2717xl0 -11 


- 

4226.997 

- 

4440.484 

cyc-C 3 Hg a) 

_ 

1.4881X10" 11 

_ 

402.4295x 10 -11 


- 

4781.938 

- 

5161.620 

C 3 Hg a) 

_ 

5.5707xl0 -13 

_ 

597.9850xl0 -13 

- 

3537.025 

- 

7118.782 

iso-C 4 H^g ) 

_ 

2.7970xl0 -14 

_ 

208.3210xl0" 14 



1598.004 


7103.169 


(a) When the potential energy surfaces are positive, the Langmuir constants are negligible and the 
parameters Ax\t and Bx,t are not given. 


experimental ones. 

To extract a quantitative value for our comparison, we use the following least square deviation pa¬ 
rameter a defined as: 


a = 


\ 


N 

(Ccal - Cy ) 2 
i —1 _ 

TV- 1 


(9) 


where C ca \ are the values calculated at different temperatures in this work and Cy are the values obtained 
experimentally [4] or from the spherical model as often used in the literature. TV is the number of 
temperatures considered (TV = 10 in Table 5). We find values of a equal to 1.6 xl0~ 2 Pa -1 and 3.3 
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Table 5: Calculated and experimental Langmuir constants (Pa 4 ) for ethane trapped in large cage of 
structure si. Experimental values are obtained from reference [3]. 


T(K) 

C ca i. this work 

C e xp. 

C ca i. Spherical model 

200 

5.3 x 10 " 2 

9.9 x 10 " 2 

14.7 xlO " 2 

210 

2.0 x 10" 2 

3.2 x 10 " 2 

5.2 xlO " 2 

220 

7.8 x 10 " 3 

12.0 xlO " 3 

20.2 xlO " 3 

230 

3.4 x 10 " 3 

4.6 xlO " 3 

8.5 xlO " 3 

240 

1.6 xlO " 3 

1.9 x 10 ~ 3 

3.9 xlO " 3 

250 

7.7 xKT 4 

8.4 xlO " 4 

18.6 xlO " 4 

260 

4.0 xlO " 4 

4.1 x 10 " 4 

9.5 xlO " 4 

270 

2.2 x 10" 4 

2.0 xlO " 4 

5.1 xlO " 4 

280 

1.3 xlO " 4 

1.1 x 10" 4 

2.9 xlO " 4 

290 

7.5 x 10 " 5 

6.1 xlO " 5 

16.7 xlO " 5 


Table 6: Calculated and experimental Langmuir constants (Pa 4 ) for cyclo-propane trapped in large cage 
of the si structure. Experimental values are obtained from reference [3|. 


T{ K) 

C C ai. this work 

n 

'-''exp. 

C C ai. Spherical model 

240 

6.7 xlO " 3 

9.5 xlO " 3 

0.12 xlO " 3 

250 

3.0 xlO " 3 

3.8 xlO " 3 

0.06 xlO " 3 

260 

1.5 xlO " 3 

1.7 xlO " 3 

0.04 xlO " 3 

270 

7.3 xlO " 4 

7.6 xlO " 4 

0.23 xlO " 4 

280 

3.9 xlO " 4 

3.8 xlO " 4 

0.15 xlO " 4 

290 

2.2 xlO " 4 

2.0 xlO " 4 

0.10 xlO " 4 
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xlO -2 Pa^ 1 , respectively. Moreover, note that in our calculations, the minimum of the potential energy 
corresponds to ethane’s center of mass located at 0.12 A from the center of the large cage of structure si, 
i.e., a value that is consistent with observations from single crystal X-ray diffraction studies that place 
the center of mass at 0.17 A from the center of the large cage flTTl . 

As for cyclo-propane (values given in Table 6 (N = 6)), we calculate values of a equal to 1.3 xlO -3 
Pa^ 1 and 3.3 xlO -3 Pa~ 3 , respectively. 

To summarize, when compared to experimental results, the values of the Langmuir constants calcu¬ 
lated here for ethane show a better agreement with experimental data than those calculated from the 
spherical model, although they are of the same order of magnitude. 

Nevertheless, for cyclo-propane Table 6 shows that values obtained from a simple spherical model 
disagree with experimental data and then they are certainly not correct. 

Therefore, the comparison of values calculated in this work with values issued from experimental work 
provides firm and sufficient grounds to validate the results obtained here for ethane and cyclo-propane. 

5 Conclusions 

In the present work, calculations of the Langmuir constants and their temperature dependence for simple- 
guest clathrate hydrates have been performed for nineteen gas species using the van der Waals and 
Platteeuw model and an all-atom approach for calculating the interactions between guest and water 
molecules in the clathrate cages. This approach accounts for the atomic character of the guest species 
and for the non-spherical water environment. 

Then, the temperature dependence of the Langmuir constants in the range 50 K - 300 K is given in 
the form of a van’t Hoff expression with parameters obtained from a fit to calculated values. 

This simple expression could thus be easily applied in various situations, especially when experimental 
data are not available. For example, it can help planetologists in the determination of the fractional 
occupancies of gas species trapped in the clathrates that are suspected to exist in various solid bodies 
of the Solar System. Indeed, this information is mandatory to analyze the atmospheric compositions 
of planets and satellites like Mars and Titan, for example, and thus to better understand their way of 
formation. Unfortunately, the temperatures usually considered in clathrate experiments are often far from 
the temperature range of interest for Planetology. As a consequence, calculations of clathrate occupancies 
is a preliminary step for any formation scenario that would consider clathrate influence. In that purpose, 
using the van’t Hoff expression given here would be much more simple than the usual approaches. 
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